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Abstract 

By simulating a uniform electric field on a lattice and measuring the change in the rest mass, we 
calculate the electric polarizability of neutral mesons and baryons using the methods of quenched 
lattice QCD. Specifically, we measure the electric polarizability coefficient from the quadratic 
response to the electric field for 10 particles: the vector mesons and K*^; the octet baryons 
n, TP, Aq, A^, and rP; and the decouplet baryons A'^, S**^, and S**^. Independent calculations 
using two fermion actions were done for consistency and comparison purposes. One calculation 
uses Wilson fermions with a lattice spacing of a = 0.10 fm. The other uses tadpole improved 
Lusher- Weiss gauge fields and clover quark action with a lattice spacing a = 0.17 fm. Our results 
for neutron electric polarizability are compared to experiment. 

PACS numbers: ll.15.Ha, 12.38.Gc, 13.40.-f 
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I. INTRODUCTION AND REVIEW 



Electric and magnetic polarizabilities characterize the rigidity of both charged and un- 
charged hadrons in external fields and are important fundamental properties of particles. 
In particular, the electric polarizability of a hadron characterizes the reaction of quarks to 
a weak external electric field and can be measured by experiment via Compton scattering. 
In this paper we describe a lattice quantum chromodynamics (QCD) calculation of neutral 
hadron electric polarizabilities using an external field method. The goal of Monte Carlo lat- 
tice QCD is to extract fundamental, measurable quantities directly from the theory without 
model assumptions. Learning about such aspects of particles tests our understanding and 
formulation of the underlying theory and makes new aspects and predictions of the the- 
ory subject to experimental verification. Such calculations can give insights on the internal 
structure of hadrons and the applicability of chiral perturbation theory to various low energy 
aspects of lattice QCD. 

Conceptually, electric and magnetic fields, even from a single photon, will distort the 
shape of a hadron, thereby affecting the internal energy and thus the mass. The electric and 
magnetic polarizabilities are defined as the coefficients of the quadratic electric and magnetic 
field terms in the mass shift formula {h = c = 1 gaussian units): 

Am = -laE'-^(3B\ (1) 

where a {(3) is the electric (magnetic) polarizability that compares to experiment. When 
the terms in Eq.(^ are viewed as the non-relativistic interaction Hamiltonian, one obtains 
the polarized on-shell Compton scattering cross section for neutral particles: 

^ = (4^? (5) I " ■ + *') ■ (''^* ' 

uji^2 and ki^2 are the initial, final photon angular frequencies and wave vectors and ei^2 are 
the polarization vectors. (The quantity 102! oj\ is a recoil factor and can be ignored nonrel- 
ativistically.) This allows hadron polarizabilities to be measured in scattering experiments. 
(For charged spinless particles one must also add the Thompson scattering amplitude from 
the particle's charge and, if virtual, a charge radius term before calculating the cross section. 
See Eq.(ll) of Q].) Compilations of experimental results for neutron electric polarizability 
have been given in . Although the polarizabilities of the other particles investigated 
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here have not been measured, we hope that the comparison of the results from the various 
types of mesons and baryons investigated will give insights on their relative rigidity and 
structure. 

The lattice calculation of electric polarizabilities began with the paper by Fiebig et al. ^ 
using the staggered fermion formulation of lattice quarks. An external electric field was 
simulated on the lattice and mass shifts measured directly for the neutron and neutral pion. 
Although the simulation errors were large, the neutron electric polarizability extracted there 
is now seen to be in remarkable agreement with recent experiments. Lattice four point 
function techniques have also been designed to extract neutral or charged particle electric 
polarizabilities [l| (chiral symmetry can sometimes be used to reduce the calculation to 
two point functions P), but these methods are more difficult to carry out on the lattice. 
Early results of the present study have been reported in ^. See Q for preliminary results 
of a companion calculation of the magnetic polarizability of both charged and uncharged 
hadrons. 



II. LATTICE DETAILS 

The clover part [7j of this calculation uses the tadpole- improved clover action with cou- 
pling constant csw = l/^o, where Uq is the average plaquette. We use the zero- loop, tadpole 
improved Liisher- Weiss gauge field action on a quenched 12^ x 24 lattice with j3 = 7.26 
(^clover = 0.17fm). In both Wilson and clover cases the gauge field was thermalized by 
10,000 sweeps (quasi-heatbath with overrelaxation) and then saved every 1000. We have 
used 100 configurations in the clover case. In the Wilson case the lattice is 24^ and we used 
109 configurations with standard gauge field action with j3 = 6.0 (awuson = 0.1 fm) j^. Our 
quark propagator time origin, t = 0, was chosen to be the third lattice time step for clover 
fermions and the second time step in the Wilson case. Point quark sources were constructed 
for the zero momentum quark propagators. The standard particle interpolation fields for the 
octet i9| (non-relativistically non-vanishing) and decouplet [l^ baryons were used. Periodic 
boundary conditions in the spatial directions and free or Dirichlet boundary conditions for 
the time links on the lattice time edges were used. 

As is usual in lattice calculations, the mass of a hadron can be calculated from the 
exponential time decay of a correlation function. By calculating the ratio of the correlation 
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function in the field to tliat witliout tlie field, we have a ratio of exponentials which decays 
at the rate of the mass difference. According to Eq.(^, the electric mass shifts should be 
parabolic, negative mass shifts giving positive polarizability coefficients and positive shifts 
giving negative coefficients. We use four values of the electric field to establish the parabola. 
Following the technique in |^, we will average the correlation functions over E and —E in 
order to remove the spurious linear term in the parabolic fits. We include the uniform static 
E- field as a phase on the gauge links in a particular direction (with fermion charge q = Qe): 

There is a discontinuity in the electric field at the decoupled lattice time boundary under this 
formulation, which should not be a problem as all of our correlation functions are measured 
far from the time edges of our lattices. Since the electric field is linearized in the continuum, 
we used the linearized form on the lattice. Fiebig et al. found no significant difference 
between the exponential and linearized formats for similar electric fields. 

On the lattice we have an exact SU(2) isospin symmetry. This symmetry is broken 
by the different electric charges on the u- and d-quarks in the presence of the external 
electric field. Consequently the vr" and / =1 particles can mix with / = glueballs and 
disconnected quark loops (self-contractions of the interpolation fields) can propagate these 
particles. Disconnected loop methods have been considerably improved in recent years 3| ; 
however, calculation of these diagrams is still extremely time consuming. In this paper, we 
will ignore the effect of the disconnected loops. 

Including the static electric field as a phase on the links affects the Wilson term, but 
not the clover loops. This is clear since the conserved vector current derived by the 
Noether procedure is identical for the clover and Wilson actions. In units of lO^^e^^a"^, 
the electric field took the values ±1.08, ±2.16, ±4.32, and ±8.64 via the parameter 
7] = o?QeE in Eq. Q for both the clover and Wilson cases. (The rj values were 
±0.00036, ±0.00072, ±0.00144, ±0.00288, and ±0.00576, the same as in Ref. Q.) In con- 
ventional units the smallest electric field is approximately 7.4 x 10^^ N/C in the clover case 
[(.17/.l)^ = 2.35 times larger in the Wilson case], which is about the same electric field 
strength .26 fm from a d-quark. Although these are huge fields, we will nevertheless see that 
the lattice mass shifts are small and that there is no evidence of E^ or higher terms in the 
electric field fits. We will in fact check that a drastic reduction in the field does not change 



4 



the polarizability coefficients. 

We did the calculation with both Wilson and clover fermions in order to test the consis- 
tency of our results as well as to set benchmarks for future polarizability calculations. We 
do not expect that the two formulations will agree with one another throughout the mass 
range investigated. As in any lattice calculation, there are many sources of systematic error, 
including quenching, finite volume, and finite size errors. The finite volume errors for the 
Wilson case {L ^ 2.4 fm) should be slightly smaller than for clover (L ~ 2 fm); however, 
discretization errors should be smaller for the clover case. In addition, to achieve physical 
results, a chiral extrapolation to physical quark masses is necessary. The quenched chiral 
regime has been estimated to extend to pion masses of only about 300 MeV for the 
pseudoscalar decay constant, fp, and the ratio m^/m, where m is the quark mass, although 
this does not exclude the range being larger for other quantities. Our pion masses are almost 
certainly too large to get into the chiral regime and we do not attempt a chiral extrapolation 
here. However, we would expect the agreement between the two calculations to improve at 
our lower pion masses, which is what is seen. Our pion masses here range from about 1 GeV 
to about 500 MeV; we achieve the smallest pion mass in the case of clover fermions, 483 
MeV. 

With = 0.1232(1) in the clover case, we used six values of k = 0.1182, 0.1194, 0.1201, 
0,1209. 0,214. and 0.1210. wh.ch .ough.v co„espo.d to ^ 200. 150. 120. 00. 70. and 

50 MeV. The critical value of k for the Wilson case is 0.157096(28)lf The values 

of K used in the Wilson calculation were 0.1515,0.1525,0.1535,0.1540,0.15454, and 0.1555, 
which correspond to approximate quark masses of rriq ~ 232, 189, 147, 126, 106, and 65 
MeV. We choose k = 0.1201 to represent clover strange quarks and k = 0.1540 to represent 
Wilson strange quarks. We used multi-mass BiCGStab as the inversion algorithm for both 
the clover and Wilson cases. The number of quark inversions per gauge field and quark 
mass was 11: 10 values of rj in Eq.Q associated with 4 nonzero electric fields for both the 
u and d-quarks, plus the zero field inversion. Our Wilson and clover calculations were done 
completely independently from one another. 
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FIG. 1: The effective lattice mass sliift, AMa, for the vr^ for the six quark masses in the clover 
study for the smallest electric field value. Error bars are shown only on the k = 0.1182,0.1219 
values. 
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FIG. 2: The effective lattice mass shift, AMa, for the 7r° for the six quark masses in the Wilson 
study for the smallest electric field value. Error bars are shown only on the k = 0.1515,0.1555 
values. 
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FIG. 5: Same as Fig. 1 (clover case) except for the neutron. 
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FIG. 6: Same as Fig. 2 (Wilson case) except for the neutron. 
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FIG. 8: Same as Fig. 2 (Wilson case) except for the A*^. 
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III. RESULTS 



We start with an examination of the effective mass shifts of the four non-strange particles 
in this study, 7r°, p", neutron, and A°, as a function of lattice time; see Figs-^IHl The effective 
mass shift for each particle is defined to be 

- SI- 

AMaW^,„(^). (5) 

where Go(^) is the particle propagator without an electric field and GE{t) is the value with 
the field. The notation AMa {t) indicates this value is being associated with the time point 
t (measured from the hadron source) in the graphs. These figures are used to guide us in 
the choice of optimal propagator time points in the fits. Just as single exponential behavior 
should emerge for each particle channel in Euclidean time, the ratio of particle propagators 
must also become single exponential. This means, from Eq.(jS)), that the effective mass shift 
plot should become time independent. As in most lattice simulations, particle propagators 
become increasing noisy as lattice time increases, so the mass shift data will eventually 
become dominated by statistical errors at large time separations. We will examine the 
results at the lowest electric field value; the effective mass shifts at larger fields we will soon 
see is simply scaled with the E'^ value. 

Figs. H El El and [7| give the non-strange hadron mass shifts for six quark masses of clover 
fermions; Figs. jSHHEl and |H1 give similar shifts for six Wilson masses. Error bars computed 
with the jackknife technique are shown on the largest and smallest quark mass results to 
give an idea of the statistical errors. We show only those lattice time points which we feel 
have meaningful Monte Carlo values in each of these graphs for clarity. Our pion result, 
Figs, m and |21 is surprising. We never do see the expected time independence of the mass 
shift for any value of the quark mass. Our results in Figs. ^ and |21 are plotted out to the 
final time step (20 in the clover case, 21 in the Wilson) to show the complete time history. 
There is apparently a dip in the time behavior for clover fermions near t = 15, but our 
error bars are too large at this point to confirm this as a plateau. The Wilson case also 
has no convincing plateau. We note that Fiebig et al. isolated a small signal for the 
neutral pion in their calculations. They used staggered fermions, which have a reminant 
exact chiral symmetry. Of course the chiral symmetry is broken for both Wilson and clover 
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fermions, and this could be the crucial difference in the calculations. It is possible that 
the neglect of the disconnected diagrams could be responsible for the bad pion behavior. 
The prediction from chiral perturbation theory is that a^o should be small and negative: 
a^o = -0.35 ± .1 X 10-^ fm^ Q. 

The other particle channels studied behave conventionally. Figs. 121 and 0] show the time 
evolution of the p° signal, which is very noisy. Although the signal decays quickly, a short 
time plateau is apparent in the data for both clover and Wilson fermions. Similar behavior 
is seen in Figs. |S] and IHl in the case of the neutron, for which the statistical signal is better. 
(Fig. ini shows the clover neutron mass shifts out to t = 16, and will be discussed below.) A 
plateau also appears in Figs. [7| and |H1 for the A°. The time plateaus appear at larger time 
steps for the Wilson data, as would be expected from the smaller value of the lattice spacing. 
The optimal time fit ranges are a compromise between statistical errors and systematic 
time evolution. Using the effective mass shift results and chi-square values as a guide, and 
assuming single exponential time fits of the ratio, R{t), we fit propagator ratio data for 
a given particle, mass, and electric field across three time steps. Although we find that 
it is often possible to fit at smaller times for larger k, we prefer to choose time plateau 
ranges independent of the quark mass being studied in order to avoid introducing spurious 
systematic effects. 

For clover data, a time plateau in the mass shift data for the neutron. Fig. begins about 
time step 5. We take this as typical of the octet baryons and fit the others in this same 
time range. The A° in Fig. [3 is noisier than the neutron but also has a plateau in the same 
region, evident for the heavier masses; we fit the same time range for the decouplets as for 
the octets. The p° in Fig. El is even noisier than the A*^, but has a plateau with acceptable 
chi-square values slightly further from the time origin. Our final choices for the optimal 
time plateaus are t = 5 — 7 for the octets and decouplet particles, and t = 6 — 8 for the p° 
and K*^. Similar considerations guide our choices of time plateaus in the Wilson case, and 
we choose to fit time steps 14 — 16 in the case of octet baryons (see Fig. ^ and 9—11 for 
the decouplets (Fig. |H1), which again are noisier. We also fit the p° (Fig. ^ and with 
time steps 13 — 15 and 10 — 12, respectively. The A^, the flavor-singlet octet member, is 
a special case. It is much noisier than the other octet members and we were not able to 
extract results from the Wilson simulations with the same time plateau as the other octet 
baryons. In the clover case, the polarizability is the largest of all the 10 particles listed, but 
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also with very large error bars. 

Figs. inK;b,c, and d are a striking demonstration of the E"^ dependence of all of our mass 
shifts. These figures show the mass shifts, Eq.(jH)), for the neutron for all quark masses and 
all four nonzero electric field values for clover fermions. Since our nonzero electric field 
values differ by factors of 2, we have rescaled the mass shifts for each of these by a factor 
of 4. Clearly, the time evolution of the mass shifts is essentially identical in these graphs 
- even the error bars scale with the factor of 4. This means that quadratic dependence on 
the electric field is assured for any choice of time fit range. The other particles, including 
the pion, display the same type of behavior, and this is seen for both clover and Wilson 
fermions. When we fit both E"^ and E'^ coefficients, we find no evidence of E^ or higher 
terms and the chi-squared values on our electric field parabola fits are quite small. 

Although Fig. El leaves little doubt about the i?^ dependence in our electric field range, 
we decreased the field in the clover calculation by a factor of 10 for the first 20 configurations 
and reexamined the mass shift plots to see if the parabola shapes changed. The shifts were 
100 times smaller than before, but the time behavior of the plots was almost identical to 
the original field strength configurations for the same 20 configurations. 

Figs. ^1 and ^2 show examples of fit parabolas and error bars for the neutron for all six 
values of quark mass. The time fit range of the propagators is 5 — 7 in the clover case and 
14 — 16 in the Wilson case. Again, Fig. 121 guarantees a quadratic field dependence. The 
averaging procedure described earlier over E and —E has removed odd terms in the electric 
field, and we just saw that there is no evidence of E^ or higher terms in the results. Figs. IT^ 
EOlshow the pion mass dependence of the extracted values of the electric polarizability, a, for 
both clover and Wilson fermions. We do not attempt to extract the 7r° or polarizabilities 
because of the problems encountered in finding a fitting plateau. Our clover results for the 
other non-strange particles are very similar to the preliminary results given in Q]. We will 
discuss the particles by groups: first the mesons, next the octet baryons, and finally the 
decouplet baryons. All are plotted on the same vertical scale (10~^ fm'^) so that error bars 
may be compared. Tabled tabulates the final results for the electric polarizability coefficient 
for various hadrons for clover fermions; Table ITTl gives the same results for the Wilson case. 

The results for the p° meson are given in Fig. ^1 There are definite incompatibilities 
in the Wilson and clover signals at the larger pion masses. However, as discussed in the 
last section, we do not necessarily expect the Wilson and clover formulations to agree at 
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FIG. 9: The neutron effective lattice mass sliifts, AMa, evaluated for 4 values of the external 
electric field, Ei = —1.08 x W~^e~^a^'^, -2Ei, 4Ei, and -8Ei in (a)-(d) respectively, for six quark 
mass values. Error bars shown only on the k = 0.1182,0.1219 values. The vertical axis scale 
increases by a factor of 4 in each case, resulting in strikingly similar figures. 
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FIG. 10: The effective lattice mass shift, AM (in units of a^il^^^) as a function of electric field in 
units of 10~^e~^a^^yj,j. for the neutron for the six values of quark mass for clover fermions. The 
time fit range of propagators is 5-7. 
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FIG. 11: The same as Fig. 10, but for Wilson fermions. The time fit range of propagators is 14-16. 
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the larger pion masses. At the smallest pion masses, the results seem to be becoming more 
compatible, although the error bars in both cases are getting quite large. The result at 
our smallest quark mass is in the range ~ 5 — 10 x 10~^ fm^. The results for the K*^ are 
given in Fig. E| In both calculations there is a significant reduction in the polarizability 
coefficient for K*^ as compared to the p°; the reduction is slightly larger in the Wilson case. 
The results for the K*^ do not seem to be converging as well as for the at the smallest 
quark mass. Note that the fit in the Wilson K*^ case has been moved 4 time steps closer to 
the propagator origin to achieve a signal; this could be contributing to the larger reduction 
of the Wilson results . The result here is in the range ~ 1 — 4 x 10^^ fm^. 

The octet baryons (n, A^, S°) are represented by Figs. (T^lfTTI The signals here are 
the best of all the particles for both Wilson and clover fermions. The behavior of the all 
four octet members are quite similar to one another when the Wilson or clover results are 
compared among themselves. The clover results tend to be a bit larger than the Wilson ones 
in all cases. All tend to move toward larger values as the pion mass is decreased. There is 
a convergence of the results for smaller masses. (The in Fig. El is a possible exception.) 
As we will discuss shortly, the trend of the neutron results are compatible with experiment. 
All particles in Figs. IT^HTTI end up with values in the range of ~ 10 — 15 x 10~^ fm^ at our 
smallest pion mass. We do not present a graphical result for the A°, but the values obtained 
in the clover case are given in Table HI 

The decouplet baryons (A°, are represented by Figs. I18ll20l Here the relative 

signals are not as strong as for the octets (most evident at smaller pion mass) and the 
values themselves are smaller. Like the octet case, the Wilson or clover decouplet results 
are quite similar to one another when compared among themselves. While we still have the 
usual incompatibilities at large pion masses, there is a trend for the results to approach one 
another for the smaller masses. (The largest disagreement is the S*° in Fig. OUj). The results 
at the smallest pion masses are all on the order of ~ 2 — 5 x 10^^ fm^. This is quite reduced 
from the octets, whose values range ~ 10 — 15 x 10~^ fm^ for small mass. As in the case 
of the meson, it was necessary to move the time steps toward the propagator origin 
in order to achieve a signal in the Wilson case; this could contribute to the relatively large 
disagreements here between Wilson and clover as opposed to the octet case. 

The only experimental result available to us for comparison is the neutron. Because of the 
absence of free neutron targets, actual Compton scattering experiments must use neutrons 
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FIG. 12: The p° meson electric polarizability coefEcient, a, in units 10 ^ fm^ as a function of the 
lattice pion mass squared in GeV^. 
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FIG. 13: The same as for Fig. 12, except for the K*^. 
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FIG. 14: The same as for Fig. 12, except for the neutron. The world average is shown at the 
physical pion mass squared. 
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FIG. 16: The same as for Fig. 12, except for the A^. 
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FIG. 17: The same as for Fig. 12, except for the 
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FIG. 18: The same as for Fig. 12, except for the A°. 
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FIG. 19: The same as for Fig. 12, except for the 11*°. 
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FIG. 20: The same as for Fig. 12, except for the S**^. 



bound in the deuteron or other nuclei. The most recent experiments have used either quasi- 
free deuteron Compton scattering (7 c? — > 'ypn) or elastic deuteron scatteri ng (/ yd — > 'yd). 
The Particle Data Group result from 2004 for the neutron is a„ = 11.6^2:3 UM- The result 
of a calculation using heavy baryon chiral perturbation theory for the neutron has been 
given as 13.4 ± 1.5 jl6|, where the error bars are associated with numerical integration and 
higher order contributions. Our Wilson and clover results are consistent with these values 
within Monte Carlo errors at our lowest pion masses, which is very encouraging. However, 
more work needs to be done, including chiral extrapolations of the present data, before more 
precise comparisons can be done. 

Since the polarizability coefficient scales like length cubed, it is extremely sensitive to 
the assigned lattice scale, a. Clearly, we must get this right if we are to extract experimen- 
tally meaningful values for the polarizabilty coefficient. The fact that both our clover and 
Wilson results for the neutron are tending toward the experimental result at the smallest 
pion masses is a strong indication of the correctness of our methods. Turning this around, 
this is a excellent place to set the lattice scale from experiment, once the neutron electric 
polarizability is determined with greater accuracy. 
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IV. CONCLUSIONS 



In this study we have used the methods of lattice QCD to evaluate the electric polariz- 
ability coefficient for neutral hadrons with both Wilson and clover fermions using improved 
Liisher- Weiss gauge fields. We were able to extract and compare electric polarizability co- 
efficients for the vector mesons and K*^; the octet baryons n, S", A^, and and the 
decouplet baryons A*^, and S*". The was extracted only in the clover case. As a rule, 
the Wilson polarizability results are smaller than the clover; however, the two calculations 
tend to converge at the smallest pion masses studied. We find that the polarizabilities of 
the octet baryons and of the decouplet baryons behave as a group; all have similar mass 
dependencies and values when the clover or Wilson results are compared among themselves. 
This is in stark contrast to the similarity in the charge radii of the proton and charged 
delta seen in Ref. Q|. Apparently, the electric polarizability is not correlated to the overall 
electromagnetic size of the hadron. In addition, Ref. sees a nontrivial squared charge 
radius behavior for the neutral baryons, the neutron having a negative charge radius and 
the others being zero or positive, which also does not seem to be reflected in the polar- 
izability. A simple harmonic oscillator model ^ for charges q and -q in an electric field 
gives a = 2q^ / {muj'^) (more generally, one obtains a sum over the squared charges of the 
constituents), where mcu^ = is the spring constant. Assuming the constitutent mass, m, of 
the quarks is a constant, a larger energy frequency uj would explain the smaller polarizability 
of the decouplets in this model. 

The polarizabilities of the octet baryons show the best agreement between clover and 
Wilson. Outside of the A°, the polarizabilities mostly tend upward as the pion mass is 
decreased. They all have values in the neighborhood of ~ 10 — 15 x 10~^ fm^ at the smallest 
pion mass. Agreement of our results with experiment for the neutron is an encouraging sign. 
There is more disagreement in the values of the decouplet baryons. Both fermion formula- 
tions agree that the values are decreased relative to the octets, but like the K*^ relative to 
the the reduction is larger for the Wilson case than the clover case. Nevertheless, the 
clover results are tending downward toward the Wilson values, and both end up with values 
in the range 2 — 5 x 10~^ fm^ at the lowest pion mass. The greatest incompatibility in the 
calculations appears in the meson sector, especially the K**^ vector meson. A possible source 
of the disagreement between Wilson and clover vector mesons and decouplet baryons seems 
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to be the shorter time interval that was necessary to achieve a signal for these particles in 
the Wilson case. We have also seen that the pseudoscalar mesons do not have an identifiable 
mass shift plateau region for either lattice fermion formulation. 

Much theoretical work remains to be done in the field of hadron polarizability. Besides 
extending the present calculations to the chiral regime, magnetic polarizabilities can be 
measured using external field techniques for both charged and neutral hadrons. Generalized 
polarizabilities, measured in polarized photon, polarized proton scattering, are beginning to 
be studied and measured |l8|, and are candidates for lattice calculations. There is also a 
need to extend chiral perturbation theory calculations of Compton scattering to isolate the 
polarizability coefficients for mesons and the other octet and decouplet baryons so that chiral 
extrapolations of lattice data may be done. In addition, the disconnected diagram contribu- 
tion to the meson polarizabilities and should be investigated, especially in the pseudoscalar 
channel. 
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TABLE I: The electric polarizabilities from the calculation with clover action using six k values. 
The units of the electric polarizability are 10""^ fm^. The pion masses were fit on time steps 11 — 13 
from the propagator origin and are given in GeV. 



K n 1 1 89 

rxj U.-L±0^ 


n 1 1 Q/1 


n 1 9(11 


n 1 909 


n 1 91 4 


n 1 91 9 






n 814 + 007 

W.O±t: _I_ .\J\J I 


f) yqr 1 nrio 

L/. 1 tJfJ -1— .\J\J<J 


f) «qc: _L nno 


f) c:«/i _|_ nnq 

\J.iJ\J'-t _l_ .WL/i/ 






Mesons 












nt range 


«0 11/11 n 
p ii.4 ± . / 


iz.U dz i.U 


1 O O _L 1 A 

lZ.Z dz 1.4 


11. o dz Z.U 


lU. / zb o.U 


O.y dz o.D 


D to O 


K*° 4.7 ± .4 


4.8 =b .5 


4.9 zb .5 


4.8 zb .7 


4.8 zb .8 


4.1 zb 1.0 


6 to 8 


Baryon octet 














n 12.8 ± .4 


13.6 lb .6 


t A f\ \ CI 

14.0 ± .8 


14.5 ± 1.2 


14.6 ± 1.6 


14.0 ± 2.5 


5 to 7 


2j 11.9 ± .0 


ion _L T 

Iz.U ± .7 


ion n 
12. ± .9 


1 1 n 1 o 
11.9 ± 1.2 


1 1 O _L 1 

11.8 ± 1.6 


lie _L O /I 

11.5 ± 2.4 


5 to 7 


A° 12.5 ± .4 


13.0 ± .6 


13.3 ± .8 


13.6 ± 1.1 


13.9 ± 1.4 


14.1 ± 1.7 


5 to 7 


AO 30 ± 1 


32 ± 21 


40 ± 24 


52 ± 35 


55 ± 43 


53 ± 86 


5 to 7 


13. 1 ± .5 


13.8 ± .7 


14.0 ± .8 


11.2 ± 1.0 


11.5 ± 1.1 


15.1 ± 1.3 


5 to 7 


Baryon decuplet 














A" 8.8 ± .4 


8.4 ± .6 


8.1 ± .9 


7.6 ± 1.4 


6.9 ± 2.1 


5.6 ± 3.3 


5 to 7 


i;*o 8.0 ± .4 


7.2 ± .6 


6.5 ± .7 


5.5 ± 1.0 


4.5 ± 1.4 


3.5 ± 2.0 


5 to 7 


2*0 7.3 ± .5 


6.8 ± .6 


6.5 ± .7 


6.1 ± .9 


5.6 ± 1.2 


5.2 ± 1.5 


5 to 7 



24 



TABLE II: The electric polarizabilities from the calculation with Wilson action using six k values. 
The units of the electric polarizability are 10""^ fm^. The pion masses were fit on time steps 11 — 13 
from the propagator origin and are given in GeV. 



K c\ ^^^^ 

rxj \J.i-Oi-0 




n 1 '\'v\ 

\J . i-OOO 


VJ. l-O^U 


n 1 '^d'^ 
\j . ±0^0 


\J.i-000 










f) 791 _L nnc 
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fit range 


0^ 5 ± 5 


4.8 ± 0.7 


4.4 lb 0.9 


4.0 lb 1.2 
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2.8 lb 4.8 
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jy-*0 1 9 _|_ ri 9 
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1 o _i_ n o 

± . KJ _l_ \J.fJ 








1 D-l 9 
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f. 7 q _|_ fi 
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in S -1- 1 Q 






SO 7.7 ± 0.7 


8.1 ± 0.9 


8.8 ± 1.3 


9.3 ± 1.5 


10.0 ± 1.9 


11.5 =b 4.0 


14-16 


AO 8.2 ± 0.7 


8.8 ± 0.8 


9.7 ± 1.1 


10.2 ± 1.2 


10.8 lb 1.6 


13.2 lb 3.2 


14-16 


SO 9.3 ± 0.9 


9.7 ± 1.0 


10.0 ± 1.2 


10.2 ± 1.4 


10.3 lb 1.6 


10.1 lb 2.3 


14-16 


Baryon decuplet 














AO 1.7 ± 0.1 


1.6 ± 0.2 


1.5 ± 0.3 


1.5 lb 0.4 


1.6 ± 0.6 


2.3 lb 1.2 


9-11 


S*o 1.7 ± 0.2 


1.6 ± 0.2 


1.5 ± 0.4 


1.5 ± 0.4 


1.6 ± 0.5 


2.0 zb 0.9 


9-11 


S*o 1.7 ± 0.2 


1.6 ± 0.3 


1.5 ± 0.4 


1.5 lb 0.4 


1.5 lb 0.5 


1.6 lb 0.8 


9-11 
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